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Abstract 


Consider  the  lifelengths  T\,...,Tk  of  k  components  subjected  to  a  randomly 
varying  environment.  They  are  dependent  on  each  other  because  of  their  com¬ 
mon  dependence  on  the  environment.  The  parameters  of  the  model  are  the 
distribution  of  the  random  process  which  describes  the  environment  and  a  set 
of  rate  functions  which  determine  the  probability  law  of  T\, . . . ,  Tk  as  a  function 
of  the  distribution  of  the  environment.  We  find  conditions  on  the  parameters  of 
the  model  which  imply  that  7\, . . . ,7\  are  associated.  Other  conditions  which 
imply  that  Ti, ...  ,7k  have  the  multivariate  aging  properties  IHR  (increasing 
hazard  rate)  and  NBU  (new  better  than  used)  are  also  described.  Also  two  such 
models  axe  compared.  In  particular,  we  characterize  the  parameters  of  these 
models  so  that  stochastic  ordering  between  the  two  vectors  of  resulting  lifetimes 
can  be  obtained. 
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1.  INTRODUCTION 

Consider  the  lifelengths  Ti,...,Tjc  of  k  components  subjected  to  a  randomly 
varying  environment.  They  are  dependent  on  each  other  hecause  of  their  com¬ 
mon  dependence  on  the  environment.  In  the  model  introduced  by  Cinlar  and 
Ozekici  (1987)  to  handle  such  dependence,  the  cumulative  hazard  functions  of 
the  components  are  made  functionals  of  the  environment  process  and  jointly 
satisfy  a  differential  equation.  Therefore,  the  joint  probability  law  of  the  life- 
lengths  is  specified  by  the  probability  law  of  the  environment  process  X  and 
the  intrinsic  aging  rates  r,(i,  a\ , . . . ,  a^),  t  G  {1, . . . ,  k},  where  the  latter  stands 
for  the  instantaneous  failure  rate  of  the  component  t  at  an  instant  when  the 
environmental  state  is  i  and  the  intrinsic  ages  (the  cumulative  hazards)  of  the 
components  l,...,fc  are  ai,---  ,  a*  respectively.  We  shall  make  these  precise 
shortly,  in  Section  2. 

Our  aim  is  to  explore  the  dependence  of  the  lifelengths  on  the  function 
r  =  (rj,...,rfc)  and  the  process  X.  In  Section  3,  we  examine  the  effects  of 
replacing  r  and  X  by  another  function  r  and  another  process  X,  in  both  cases 
seeking  results  on  stochastic  dominance.  Also  in  that  section  is  a  characteriza¬ 
tion  of  “association”  (in  the  sense  of  Esary,  Proschan  and  Walkup  (1967))  for 
the  lifelengths  in  terms  of  the  association  of  the  process  X. 

In  Section  4,  we  consider  multivariate  aging  properties  of  the  lifelengths  con¬ 
ditioned  upon  the  history  7t  of  the  environment  until  t,  and  also,  conditioned 
upon  the  history  Qt  of  the  environment  and  failures  during  (0,  t].  In  particu¬ 
lar,  we  obtain  conditions  for  the  lifelengths  to  have  the  “multivariate  increasing 
hazard  rate”  property  with  respect  to  the  filtration  (J<)  or  (£t),  and  also  the 
“multivariate  new  better  than  used”  property,  again  with  respect  to  (7<)  or  [Qt)- 


2.  PRELIMINARIES 


In  this  section  we  give  an  overview  of  the  model  introduced  by  Qinlar  and 
Ozekici  (1987).  Throughout  here  and  the  paper,  (fi ,M,P)  is  a  complete  proba¬ 
bility  space.  We  write  R+  for  [0,  oc),  call  a  number  or  vector  a  positive  [negative) 
if  a  >  0[a  <  0),  and  call  a  function  /  increasing  [decreasing]  if  /(x  )  <  f(y)  for 
a:  <  >  y]. 

Let  ( E ,  £)  be  a  measurable  space.  Elements  of  E  are  called  the  environmental 
states.  We  suppose  that,  for  each  x  €  E,  the  singleton  {x}  belongs  to  £ .  There 
is  a  distinguished  point  in  E ,  denoted  by  6,  which  stands  for  the  state  that 
causes  no  aging.  We  let  A'  =  {A(t);t  G  R+}  be  a  stochastic  process  with  state 
space  (E,£);  it  represents  the  environment  process. 

The  set  of  all  components  is  epresented  by  K  —  k}.  We  let  A  — 

{A(t);<  6  R+}  be  an  increasing  continuous  process  taking  values  in  J2+;  its  tth 
component,  namely  A,  =  {yl,(<);t  6  is  called  the  intrinsic  age  process  of 

component  i,  it  plays  the  role  of  a  random  cumulative  hazard  function. 

We  let  S\ , . . . ,  Sk  be  independent  of  X  and  of  each  other  and  have  the  standard 
exponential  distribution  (with  mean  1).  The  lifelength  of  component  i  is  modeled 
by 

(2.1)  Ti  =  inf {i  :  At{t)  >  5,},  i  G  K, 

that  is,  the  component  i  fails  when  its  intrinsic  age  runs  over  its  “intrinsic 
lifelength”  S,.  We  write  S  —  (Si, . . . , Sk)  and  T  =  (Tj, . . .  ,Tk)  for  the  vectors 
of  intrinsic  lifelengths  and  lifelengths. 

In  this  formulation,  the  dependencies  between  the  lifelengths  and  their  joint 
dependence  on  the  environment  are  reflected  via  the  intrinsic  age  process  A. 
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Regarding  the  latter,  the  main  assumption  of  C'mlar  and  Ozekici  (1987)  is  as 
follows. 

(2.2)  HYPOTHESES,  (i)  For  each  component  i  there  exists  a  positive  measur¬ 
able  function  on  E  x  Rk+  such  that 


(2.3)  dAi{t)  =  r{(X(t),A(t))dt,  t  >  0,  i  e  K. 

(ii)  We  have  r,(x,  a)  >  0  for  each  i  G  K,  a  6  R+,  and  all  x  €  E  except  r  —  6. 
For  x  =  6 ,  we  have  rt(x,a)  =  0  for  all  t  and  a. 


The  basic  hypothesis  is  the  first  one:  the  intrinsic  age  process  A  is  a  functional 
of  the  environment  process  X.  The  second  hypothesis  is  a  regularity  condition,  it 
is  meant  to  ensure  that  (2.3)  has  a  unique  solutic  A  for  each  starting  condition; 
in  particular,  it  singles  out  6  as  the  only  state  that  causes  no  aging. 


Note  that  A  is  determined  by  X  and,  hence,  is  independent  of  the  vector 
S.  Therefore,  it  follows  from  (2.1)  and  the  independence  of  the  exponential 
variables  Si  from  each  other  that 


(2.4) 


P{T  >  t|X}  =  exp 


E 

iGK 


t  <e  id. 


This  justifies  the  term  “random  cumulative  hazard  function”  for  each  A,.  How¬ 
ever,  we  prefer  to  call  A,  the  intrinsic  age  process  of  i.  Then  (2.1)  can  be  read  as 
follows:  each  component  is  endowed  with  an  intrinsic  lifelength,  the  component 
ages  in  response  to  the  environmental  factors  in  a  manner  intrinsic  to  its  own 
function  and  nature,  it  fails  when  its  intrinsic  age  runs  over  its  allotted  intrinsic 
lifelength.  In  particular,  (2.3)  defines  r,(x,ai, . . . ,  a*)  to  be  the  intrinsic  aging 
rate  of  the  component  i  at  a  time  when  the  environment  is  in  state  x  and  the 
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intrinsic  ages  of  the  components  1  are  ai,...,ak  respectively.  It  follows 

from  (2.3)  and  (2.4)  that  we  also  have 

(2.5)  r<(x;ai,...,a*) 

=  lim  -P{T{  <  t  +  u\Tt  >  t.  X(t)  =  x,Ax{t)  =  a1,...,Ak(t)  =  ak }, 

tt|0  It 

that  is,  r,  (x,a)  is  the  hazard  rate  for  component  j  as  a  function  of  the  environ¬ 
mental  state  x  and  the  intrinsic  age  vector  a.  We  write  r  =  (rl5. . .  ,rk)  and  call 
it  the  intrinsic  aging  rate  function.  Aside  from  the  probability  law  of  X,  it  is 
the  only  parameter  in  the  model. 

If  r(x,a)  is  free  of  x,  then  A  becomes  deterministic  and  (2.1)  shows  that 
Tj , . . .  ,Tk  are  independent.  If  r(x,  a)  is  free  of  a,  then  A  becomes  a  fc-dimensional 
additive  functional  of  X. 

A  function  r  from  ExR1^  into  R+  will  be  called  an  intrinsic  aging  rate  function 
if  it  satisfies  Hypothesis  (ii)  of  (2.2).  Given  such  a  function  r  and  the  process 
X,  the  differential  equation  (2.3)  together  with 

(2.6)  A(0)  =  0 

specifies  the  aging  vector  A{t)  for  all  t  >  0,  and  the  latter  specifies  the  lifelength 
vector  T  via  (2.1)  from  intrinsic  lifelength  vector  S  of  standard  exponentials. 
Thus,  there  exists  a  functional  L  such  that 

(2.7)  T  =  L(X,r,S). 

The  functional  L  is  defined  implicitly  via  (2.1)  and  (2.3);  it  is  called  the  lifelength 
functional.  This  paper  is  a  study  of  the  dependence  of  L  on  its  arguments  X 
and  r. 
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3.  DEPENDENCE  ON  ENVIRONMENT  AND  AGING  RATES 

In  this  section  we  discuss  the  dependence  of  the  lifelength  vector  T  = 
L(X,r,S )  on  the  environment  process  X  and  the  intrinsic  aging  rate  function  r. 
Here,  and  for  the  remainder  of  the  section,  we  assume  that  the  state  space  E  is 
a  complete  separable  metric  space. 

(3.1)  THEOREM.  Let  r  and  f  be  intrinsic  aging  rate  functions  and  let  T  = 
L{X,r,S )  and  f  =  L(X,f,S).  Assume  that  t  —*  X (t)  is  piecewise  continuous 
and  that  either  r  or  f  is  continuous  on  E  X  Rk.  Suppose  that,  for  each  t  €  K 
and  x  €  E, 

(3.2)  a,  a  €  Rk+,a  >  a,a{  =  at  =>  r,(x,a)  >  r,(x,  a). 

Then,  T  <  f. 

(3.3)  REMARK.  Suppose  that  r  >  f  and  that  either  a  —>  r(x,  a)  or  a 
f(x,a)  is  increasing  for  every  x.  Then,  the  condition  (3.2)  of  the  preceding 
theorem  is  satisfied  and  T  <  T.  For  instance,  if  r  >  f  and  a  — ►  r(x.a)  is 
increasing  for  each  x,  then  r(x,a)  >  r(x,a)  >  f(x,a)  for  all  a  >  a,  zmd  hence 
(3.2)  holds. 

(3.4)  REMARK.  Suppose  that  r  and  f  are  as  in  the  preceding  theorem  and 
(3.2)  is  satisfied.  Suppose  that  T  =  L(X,r,S)  and  T  =  L(X,r,S)  where  X  and 
X  have  the  same  probability  law,  and  so  do  S  and  S,  and  S  is  independent 
of  X  (as  S  is  of  X).  Then,  the  conclusion  of  the  preceding  theorem  is  that  T 
is  dominated  by  T  stochastically,  that  is,  Ef(T)  <  Ef(T)  for  every  increasing 
function  /  from  R+  into  R+ . 

Proof  of  Theorem  (3.1).  Fix  r  and  r.  Suppose  first  that  (3.2)  holds  with  a 
strict  inequality:  r\(x,  a)  >  fi(x,a).  Let  A  be  the  solution  of  (2.3)  and  let  A  be 


G 

the  solution  of  (2.3)  with  f  replacing  r,  both  with  .4(0)  =  .4(0)  =  0.  Now,  T  is 
defined  by  (2.1),  and  f  is  defined  by  (2.1)  with  A  replacing  A.  Thus,  to  show 
that  T  <  T,  it  is  sufficient  to  show  that  A  >  A.  Or,  equivalently,  it  is  sufficient 
to  show  that  the  random  variable 

(3.5)  r  =  inf (t  :  Ai(t)  <  A,(<)  for  some  *} 

is  equal  to  +oo  identically. 

Since  >1(0)  =  A(O)  =  0,  we  have  r  >  0.  Suppose  for  the  moment  that  r(u>)  =  t 
(where  t  <  oc)  for  some  outcome  u  £  fi.  Fix  that  w  and  simplify  the  notation 
by  putting 

(3.G)  x  —  A\(cj),a  =  A(oj,t),a  -  A(iu,t). 


In  view  of  (2.3),  the  processes  A  and  A  are  continuous,  and  the  assumed 
finiteness  of  t  =  r(w)  implies  the  existence  of  i  £  K  and  of  a  decreasing  sequence 
(t„)  C  R+  with  limit  t  such  that 

(3.7)  ^4, (u>, <„)  <  Ai(u,tn)  for  all  n. 

Moreover,  by  the  continuity  of  A  and  A,  we  must  have 

(3.8)  a  >  a,  a*  =  d< 


The  differentiability  of  Ai  and  A,  ensured  by  (2  3)  implies  that,  since  a,  =  a ,, 

1 


lim 


■[A,(w,tn)  A»(w,tn)j 


n— >oo  tn  —  t 

=  lim  -(Aj(u;,tn)  -  Ojj  -  lim  ——  - 

n->oo  tn  —  t  n— >oo  tn  —  t 

=r,(  x,a)  —  f<(x,d). 


(3.9) 


|A,(w,tn)  ujj 


In  view  of  (3.8)  and  the  assumed  strictness  in  condition  (3.2),  this  is  strictly 
positive.  But,  the  first  member  of  (3.9)  must  be  negative  in  view  of  (3.7).  This 
contradiction  shows  that  -(:*,’)  cannot  be  finite. 


Now  relax  the  assumption  of  strictness  in  (3.2)  but  assume  for  a  moment 
that  t  — ♦  X(t)  is  continuous.  Then  t  — ♦  r(X(t),A(t))  or  t  -*  f  (X(t),  A(f )) 
is  continuous  by  the  hypothesis  that  either  r  or  r  is  continuous.  Suppose  the 
former.  Fix  an  e  >  0  in  jO,  oc)k  and  define  ~  r  +  e/n,n  =  1,2,...  .  Let 
be  the  solution  of  (2.3)  with  replacing  r  and  with  Aln)(0)  =  0.  Then  by 
the  previous  argument  >  A.  By  the  continuity  of  r,  A(n)  — ►  A.  Therefore 
A  >  /I  in  this  case. 

The  proof  for  the  case,  in  which  t  — ►  f(X(f),  A[t))  jrathcr  than  t  — ► 
r(x{t),  A[t)\  is  continuous,  is  similar.  The  above  argument  can  be  used  except 
that  r  is  replaced  by  f =  f(l  -  ejn)  where  e  <  1  is  fixed.  This  definition  of 
ensures  that  r^n)  >  0  as  required  in  (2.2)  (ii). 

If  t  — >  X{t)  is  piecewise  continuous,  then  let  tL,t2.  ■  ■  ■  be  the  successive  iump 
times  of  X.  On  each  interval  (^,t£+i),  X  is  continuous  and  the  previous  argu¬ 
ment  can  be  applied  to  each  such  interval  to  show  that  A  >  A.  Ij 

Dependence  on  environment 

For  the  remainder  of  this  section,  we  assume  that  the  state  space  E  is  a 
partially  ordered  Polish  space  (a  complete  separable  metric  space  with  a  closed 
partial  ordering) .  Then,  the  space  D  =  D(R+ ,  E )  of  right-continuous  left-limited 
functions  fromi?^  into  E  is  again  a  partially  ordered  Polish  space.  A  functional 
g  :  D  R+  is  said  to  be  increasing  if  w  <  w  implies  p(ui)  <  <?(u))  for  all  paths 
w,ii>  €  D,  where  <  denotes  the  partial  ordering. 


Let  X  and  X  be  processes  with  paths  in  D.  Then,  A'  is  said  to  dominate  X 
stochastically  provided  that 

(3.10)  Eg( X)  >  Eg( X) 

for  every  Borel  measurable  increasing  functional  g  on  D.  This  is  obviously  the 
case  if  X(u.\t)  >  A for  all  u  €  fl  and  t  £  R+.  More  generally,  if  X 
stochastically  deminates  A',  then  it  follows  from  Theorem  1  of  Kamae,  Krengel, 
and  O’Brien  (1977)  that  A'  and  X  can  be  “put  on  the  same  probability  space 
so  that  one  dominates  the  other  path  by  path”.  More  precisely,  it  is  possible  to 
construct  a  new  probability  space  (IV,  £,Q)  and  stochastic  processes  Y  and  Y 
defined  on  (Wr.  C/,Q)  and  having  paths  in  D  such  that  Y(u,t)  >  T(w,t)  for  all 
u>  £  1L  and  t  1 11  X  and  )'  have  the  same  probability  law,  and  X  and  Y  have 
the  same  probability  law.  Of  course,  the  new  probability  space  can  be  enlarged 
to  accommodate  k  independent  standard  exponential  variables  independent  of 
Y  and  Y .  These  remarks  will  be  useful  in  simplifying  the  proof  of  the  following 
theorem,  which  reduces  to  Theorem  (3.1)  when  A"  =  X. 

Let  D  =  D[R^.,E)  be  the  set  of  functions  in  D{R+,E)  which  are  piecewise 
continuous. 

(3.11)  THEOREM.  Let  X  and  X  be  processes  with  paths  in  D,  let  r  and  r 
be  intrinsic  aging  rate  functions,  and  let  S  and  S  be  fc-vectors  of  independent 
standard  exponential  variables  independent  of  X  and  X  respectively.  Assume 
that  r  and  r  are  continuous  on  E  x  R+.  Suppose  that 

i)  X  dominates  X  stochastically, 

ii)  i  — >  r(x,  a)  is  increasing  for  every  a  (or  x  — »  f(x,  a)  is  increasing  for  every 
o),  and  the  condition  (3.2)  holds  for  every  t  £  K  and  x  €  E. 
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Then.  T  —  L[X,r,S)  is  stochastically  dominated  by  T  =  L(X,r,S). 

Proof.  In  view  of  the  foregoing  remarks,  by  moving  onto  a  new  probability 
space  if  necessary',  we  may  and  do  assume  that  S  =  5  and  X (w\  t)  >  X[u:,t)  for 
all  w  and  t. 

Let  A  be  as  before,  and  define  A  as  the  solution  of  (2.3)  with  X  and  r  replaced 
by  X  and  f,  A(0)  —  A(Q)  =  0.  As  in  the  proof  of  Theorem  (3.1),  it  is  sufficient 
to  show  that  A  >  A,  or  equivalently,  that  r  defined  by  (3.5)  is  equal  to  +  oo 
identically. 

First  assume  that  (3.2)  holds  with  a  strict  inequality'.  Let  r  be  defined  by 
(3.5)  and  suppose  again  that  t(lj)  =  t{t  <  oc)  for  some  y£!l.  Pick  i  €  K  and 
(tn)  c  so  that  ( tn )  decreases  to  t  and  (3.7)  holds.  With  the  notations  (3.6) 
supplemented  by  i  =  ,\'(u>,t),  (3.9)  becomes 

(3.12)  lira  1  [A^cu,^)  -  A,(uMn)j  =  r,(x,a)  -  f,(j,a). 

n  ►oo  t  n  t 

Since  A"  dominates  X,  we  have  x  =  Ar(w,t)  >  X(cu,t)  =  x.  Thus,  the  condition 
(ii)  implies  that 

r,(i,a)  >  Ti{x,a)  >  fj(x,a) 

if  x  — ♦  r(x,a)  is  increasing  (and  r,(x,  a)  >  f,(x,a)  >  f,(x,a)  if  x  — >  f(x,a)  is 
increasing).  It  follows  that  the  right  side  of  (3.12)  is  strictly  positive.  But  from 
(3.7)  it  is  seen  that  it  is  negative.  Hence  r(u>)  cannot  be  finite. 

The  extension  of  the  above  argument  to  the  case  in  which  strictness  in  (3.2) 
is  not  assumed  can  be  done  as  in  Theorem  (3.1).  || 

The  preceding  proof,  with  r  —  f,  yields  the  following  technical  result  regarding 
the  lifelength  functional  L. 
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(3.13)  COROLLARY.  Suppose  that  a  — >  r(x,  a)  is  increasing  for  every  x  6 
E,  that  x  — >  r(x,a)  is  increasing  (respectively,  decreasing)  for  every  a  €  J2*.  and 
that  r  is  continuous  on  ExRk_ l .  Then,  w  — ►  L(w,r,  s),  is  decreasing  (respectively, 
increasing)  in  w  G  D  for  fixed  r  and  s. 

4.  ASSOCIATION  OF  LIFELENGTHS 

Let  Z i, . . .  -Zm  be  random  variables  taking  values  in  Rn.  Then,  they  are  said 
to  be  associated  provided  that  the  vector  Z  =  (Zj, . . . ,  Zm)  satisfy 

Cov(g(Z),  k(Z))  >0 

for  all  increasing  functions  g,h  :  RmXn  — ►  R  for  which  the  covariance  exists. 
A  stochastic  process  Z  -  {Z{t)\t  €  JR+}  with  state  space  Rn  is  said  to  be 
associated  in  ti).  j  if  Z(t  i), . . . ,  Z(tm)  are  associated  for  all  integers  m  >  1 
and  times  ti,...,tm  G  R+.  Our  aim  in  this  section  is  to  show  that,  if  the 
environment  process  X  is  associated  in  time  and  certain  conditions  hold  for  the 
aging  rate  function  r,  then  the  lifelengths  Ti,...,Tjt  are  associated.  We  refer 
to  Esary,  Proschan,  and  Walkup  (1967),  Barlow  and  Proschan  (1975),  Arjas 
and  Norros  (1984),  Shaked  and  Shanthikumar  (1987)  and  references  therein  for 
the  usefulness  of  the  concept  of  association  for  lifelengths,  and  to  Barlow  and 
Proschan  (1976)  and  Harris  (1977)  for  examples  of  processes  associated  in  time. 

For  the  purposes  of  this  section  we  assume  that  the  environment  process  X 
takes  values  in  E  =  Rn  and  its  paths  belong  to  D  as  in  the  preceding  section. 

(4.1)  THEOREM.  Suppose  that  X  is  associated  in  time.  If  x  — *  r(x,  a)  is 
increasing  for  every  a  G  R+  (or  decreasing  for  every  a  G  R+),  a  — ♦  r(x,a)  is 
increasing  for  every  x  G  E  =  Rn,  and  r  is  continuous  on  E  x  Rk+,  then  the 
lifelengths  Tu...,Tk  are  associated. 
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Proof.  Fix  r,  suppose  that  r(x,a)  is  increasing  in  both  x  and  a.  Then,  by 
Corollary  (3.13),  the  mapping  w  — >  L(w,r,s)  from  D  into  72+  is  decreasing. 
Thus,  if  g  and  h  are  increasing  functions  from  72l  into  72,  then  —  g  o  L[w,r,  s) 
and  —h  o  L(w,r,s)  are  increasing  functions  of  w  €  D  and  we  have 

(4.2)  EgoL{X,r,s)hoL{X,r,s)  >  Eg  o  L{X,r,s)Eh  o  L{X,r,s) 


by  the  assumption  that  X  is  associated  in  time.  The  same  is  true  for  the  case 
where  x  — +  r(x,  a)  is  decreasing,  by  Corollary  (3.13)  and  the  association  applied 
directly  to  g  o  L  and  ho  L. 

Let  /i  denote  the  Ar-dimensional  standard  exponential  (that  is,  the  distribution 
of  S).  By  the  independence  of  X  and  5,  the  integral  of  the  left  side  of  (4.2) 
with  respect  to  n(ds)  is  equal  to  Eg(T)h{T).  Thus,  (4.2)  gives, 

(4.3)  Eg(T)h{T)  >  j  »{ds)Eg  o  L{X,r,s)Eh  o  L{X,r,s). 

On  the  other  hand,  it  is  obvious  that  s  — *  L(w,r,s)  is  increasing,  which  implies 
that  s  — ►  Eg  o  L(X ,  r,  s)  and  s  Eh  o  L(X,  r,  s)  are  increasing.  Since  Sx,...,Sk 
are  independent,  they  are  associated.  This  in  turn  implies  that  the  right-side  of 

(4.3)  is  greater  than  or  equal  to 

J  n{ds)Eg  o  L(X,  r,  s)  J  n(ds’)Eh  o  L( X,  r,  s') 

=  Ego  L(X,  r,  S)Eh  o  L{X ,  r,  S)  =  Eg(T)Eh(T). 


This  completes  the  proof.  |] 

In  the  preceding  theorem,  the  condition  that  X  be  associated  in  time  is  satis¬ 
fied  for  processes  X  that  have  independent  positive  increments  (e.g.  increasing 
compound  Poisson  processes,  gamma  processes,  etc.).  More  generally,  in  the 
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case  of  real-valued  processes  X,  association  in  time  holds  if  X  is  stochastically 
monotone,  that  is,  if 

(4.4)  E[g{X(t))\X(0)  =  x]  <  E[g(X(t))\X(0)  =  y] 

for  x  <  y  and  g  increasing  Borel  measurable  (see  Barlow  and  Proschan  (1976) 
and  Harris  (1977)  for  this).  Thus,  the  preceding  theorem  remains  true  if  A"  is  a 
real  valued,  stochastically  monotone  Markov  process. 

5.  MULTIVARIATE  AGING  PROPERTIES 

In  this  section  0t  will  denote  an  operator  that  shifts  the  time  origin  to  t.  In 
particular, 

(5.1)  6tTi  =  max(0,Ti  -  <),  *  -  K. 

The  following  properties  were  defined  in  Arjas  (1981). 

(5.2)  DEFINITION.  Let  ()/t)  be  a  filtration.  The  lifelength  vector  T  is  said 
to  have  a  multivariate  increasing  hazard  rate  with  respect  to  (I/t)  (abbreviated 
as  (Ht)- MIHR)  if 

(5-3)  E[f(9tT)  I  Ht]  >  E[f(OuT)  |  Mu] 

for  alii  <  u  and  all  positive  increasing  Borel  functions  /  on  R+.  It  is  said  to  have 
the  multivariate  new  better  than  used  property  with  nspect  to  (!/t)  (abbreviated 
as  (tft)-MNBU)  if 

(5.4)  E\f{T)  |  Wo]  >  E[f(0tT)  |  I/t] 

for  all  t  >  0  and  all  positive  increasing  Borel  functions  /  on  R+. 
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Two  special  filtrations  of  interest  to  serve  as  (l/()  above  are  defined  by 

(5.5)  7,  =<r(/l(0),X(a)  :«<*), 

(5.6)  C,  =  j-t  V  o(I{Tt<s}  '■  s  <  t,i  €  K ), 

which  are,  respectively,  the  history  of  environment  and  age  processes  until  t  and 
the  complete  history  of  environment,  ages,  and  failures  until  t.  Note  that  we 
allow  the  initial  ages  .<4,(0)  to  be  non-zero  random  variables. 

Our  aim  in  this  section  is  to  discuss  MIHR  and  MNBU  properties  of  T  with 
respect  to  (jt)  and  (£t)  assuming  that  X  is  a  Markov  process  with  certain 
properties.  We  start  by  some  computations  in  the  Markovian  case. 

Lifelengths  in  a  Markovian  environment 

Let  X  be  a  temporally  homogeneous  Markov  process  with  state  space  E  —  /?" . 
Suppose  that  its  paths  belong  toU,  the  space  of  all  right-continuous  left-limited 
functions  fromiZ+  into  E.  Let  A  satisfy  the  differential  equation  (2.3),  but  with 
the  initial  condition  A(0)  unspecified.  It  follows,  then,  that  the  pair  (AT,  A)  is  a 
temporally  homogeneous  Markov  process  with  state  space  E  x  Rk+.  As  is  usual 
in  the  theory  of  Markov  processes,  we  will  write 

Pxa{-}  =  P{.  j  X(0)  =  x,A[0)  =  a},  x  €  E,a€Rk+, 

and  will  write  Exa  for  the  corresponding  expectation  operator.  Note  that  Pza 
does  not  put  any  conditions  on  the  vector  S  of  standard  exponentials,  except 
that  S  is  assumed  to  be  independent  of  that  is,  of  the  process  (X,  A). 

The  lifelengths  T{  are  still  defined  by  (2.1),  which  implies  that  some  of  the  7* 
can  be  0  with  a  strictly  positive  probability.  However,  if  it  is  given  that  T{  >  0, 
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its  probability  law  is  the  same  as  that  of 

(5.7)  Ui  =  inf {f  :  At(t)  -  At{0)  >  5,},  t‘  €  K , 

which  fact  follows  from  the  independence  of  S  from  A  and  the  memory lessness 
of  exponential  variables.  The  following  is  a  precise  version  of  this  circle  of  ideas. 
Here,  and  below,  for  I  C  K  and  t;  6  R+  we  define  vj  €  R+  to  be  the  vector 
whose  i-entry  is  v,  or  0  according  as  t  6  I  ^r  not. 

(5.8)  LEMMA.  Let  /  be  a  positive  Borel  function  on  R+  and  put 

(5.9)  g(x,a,I)  =  Eza[f(T)  |  5,  >  a,-  for  t  €  /  and  Si  <  a{  for  t  €  K  -  7], 
where  x  €  E,  a  G  R+,  and  I  C  K.  Then, 

(5.10)  g(x,a,I)  =Exaf(U,). 

Proof.  Under  Pxa  we  have  A(0)  =  a.  Thus,  on  {5,  <  a,}  we  have  T,  —  0 
almost  surely,  and  (5.9)  becomes 

(5.11)  g(x,a,I)  =  Eza[f{TI)\Sl>al,  i  G  /]. 

On  the  other  hand,  on  {5,  >  a,},  we  have 

Ti  =  inf{t  :  A,(t)  >  5,}  =  inf{t  :  A,(t)  -  A,( 0)  >  5,} 

where  5,  =  S,  —  a,  since  At  (0)  =  a,  under  Pxa.  By  the  independence  of  S  from 
(AT,  A),  and  since  Si  is  exponential,  Si  —  S’,  —  a*  has  the  standard  exponential 
distribution  as  its  conditional  distribution  on  {5,  >  a*}.  It  follows  that  the 
conditional  distribution  of  T/,  given  {Si  >  aj,i  €  /},  under  Pxa  coincides  with 
the  distribution  of  Uj  under  Pxa.  Hence,  the  right  sides  of  (5.10)  and  (5.11)  are 
the  same.  || 
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(5.12)  LEMMA.  Let  /  and  g  be  as  in  Lemma  (5.8).  Let  fi  be  the  standard 
exponential  distribution  on  Rk+ ,  that  is,  fi(ds)  =  exp(— si  —  •  •  •  —  sfc)dsi  •  •  •  dsk- 
Then. 

(5.13)  Exaf(T)  =  j  n(ds)g{x,ajga)=h(x,a) 

Jr\ 

where  Isa  —  {i  €  K  :  >  a,}. 

Proof.  It  is  immediate  from  Lemma  (5.8)  by  unconditioning.  || 

The  proof  of  the  next  lemma  follows  from  the  Markov  property  of  (X,A). 
Here, 

(5.14)  R{t)  =  {ieK  :  Ti  >  t}, 
is  the  set  of  components  remaining  alive  at  t. 

(5.15)  LEMMA.  Let  X  be  a  temporally  homogeneous  Markov  process.  Let 
/, <7,  and  h  be  related  by  (5.9)  and  (5.13).  Then 

(5.15)  E[f(OtT)  ,  Qt]  =  g{X{t),  A(t),  £(*)), 

(5.16)  E[f(0tT)  |  Tt]  =  h{X{t),A{t)). 


Increasing  hazard  rates 

(5.17)  THEOREM.  Let  A"  be  a  temporally  homogeneous  Markov  process 
with  state  space  E  =  Rn.  Suppose  that 

a)  r(x,  a)  increases  in  x  and  in  a  and  is  continuous, 


b)  X  is  stochastically  monotone 
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c)  the  paths  of  .Y  belong  to  D(R+,Rn)  and  arc  increasing. 

Then,  T  has  the  ( J()-MIHR  and  (£j)-MIHR  properties. 

Proof,  i)  Let  /  be  an  increasing  function  on  Rkf  and  let  g  and  h  be  defined 
by  (5.9)  and  (5.13).  To  show  that  (5.13)  holds  with  =  7t  or  Qt,t  >  0,  it 
is  sufficient  to  show  that  (5.15)  and  (5.16)  are  decreasing  in  t.  Since  X  and  A 
are  increasing  processes  (the  assertion  on  X  is  via  the  assumption  (c))  and  R 
is  decreasing,  this  amounts  to  showing  that  g  and  h  are  decreasing  in  their  first 
two  arguments  and  g  is  increasing  in  its  last  argument. 

ii)  It  is  easy  to  sec  that  g{x,a,I )  increases  as  I  increases:  if  7  c  J  then 
Uj  <  Uj  and  /([/,)  <  f{Uj). 

iii)  Fix  a  and  I.  Since  x  — ►  r(x,a)  is  increasing,  the  random  vector  Uj  is  a 
decreasing  functional  of  X  (by  (3.13)).  By  the  assumed  stochastic  monotonicity 
of  X,  this  implies  that  g(x,  a,  I )  decreases  in  x.  Further,  in  view  of  (5.13),  h{x,  a) 
decreases  in  x. 

iv)  Fix  x  and  I.  Let  a  <  a,  and  let  A  and  A  be  the  solutions  of  (2.3)  starting 
from  a  and  a  respectively.  As  before  in  Theorem  (3.1),  using  the  assumption 
that  r(x,  a)  is  increasing  in  a,  we  see  that  A  <  A.  This  implies  that 

(5.18)  jtA[t)  =  r(JY(0,A(0)  <  r(X(t),A(t))  = 

Integrating  over  (0,t]  we  see  that  A(t)  -  A(0)  <  A(t)  -  A(0).  Thus,  in  view  of 
the  definition  (5.7)  of  U,  we  have  U  >  U  where  U  corresponds  to  A  as  U  does 
to  A.  It  follows  that  f{Uj)  >  f(Uj)  and,  since  the  law  of  U  under  Px&  is  the 
same  as  that  U  under  Pld,  we  have  that  g(x,a,I)  >  g(x,a,I).  Hence,  g(x,a,I) 
decreases  in  a. 


► 
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Finally,  fix  x  and  let  a  <  a.  Consider  the  formula  (5.13)  for  h.  For  any  s  €  12*  , 
La  =  {*  :  5,  >  a,}  D  {»  :  st  >  a,}  =  I3i,  and  hence,  g(x,a,I3a)  >  g{xta,I3&).  It 
follows  from  (5.13)  that  h(x,a)  >  h(x,a),  that  is,  h(x,a)  decreases  in  a.  |j 

Note  that  the  conditions  of  Theorem  (5.17)  imply  the  conditions  of  Theorem 
(4.1)  [See  the  discussion  following  the  proof  of  Theorem  (4.1)].  This  is  not 
surprising:  Using  ideas  such  as  in  Norros  (1985)  it  can  be  shown  that  if  T  has 
the  ((Jf)-MIhR  property  then  Tx,...,Tk  are  associated. 

In  the  preceding  proof  we  had  the  assumption  that  the  paths  of  X  are  increas¬ 
ing.  For  proving  the  generally  weaker  property  MNBU,  we  may  replace  it  with 
something  weaker. 

(5.19)  THEOREM.  Let  X  be  a  temporally  homogeneous  Markov  process 
with  state  space  E  =  JRn.  Suppose  that  the  condition  (a)  and  (b)  of  Theorem 
(5.17)  hold,  and  that 

(c')  X(0)  <  X{t)  almost  surely  for  each  t  and  the  paths  of  X  belong  to 
D(R+,Rn). 

Then,  T  has  the  properties  (Jt)-MNBU  and  (^t)-MNBU. 

Proof.  Here  we  have  lo  =  §o  =  a(X(0)/l(0)).  And,  by  the  computations  of 
Lemma  (5.8), 

E{f(T)\g0)=g(X(0),A(0),K) 

So,  we  need  to  show  (by  (c'))  that 

ff(X(0),A(0),X)>ff(X(t),A(t),12(t)) 

and  that 

A(X(O),A(0))>/i(X(t),A(0). 
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But  these  follow  from  the  proof  of  Theorem  (5.17).  |! 

Theorem  (5.19)  applies  to  “new”  components  by  setting  A(0)  =  0  with  prob¬ 
ability  one. 

Note  that  (c')  holds  whenever  E  —  72”  and  P{X{ 0)  =  0}  =  1 . 
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